Annals of Mathematics, 153 (2001), 259-296 



o 
u 



Discrete orthogonal polynomial ensembles 
and the Plancherel measure 

O . By Kurt Johansson 

(N; 

^ ■ Abstract 

D ■ 

Mh . We consider discrete orthogonal polynomial ensembles which are discrete 

■^ ! analogues of the orthogonal polynomial ensembles in random matrix theory. 

I These ensembles occur in certain problems in combinatorial probability and can 

be thought of as probability measures on partitions. The Meixner ensemble is 

related to a two-dimensional directed growth model, and the Charlier ensemble 

is related to the lengths of weakly increasing subsequences in random words. 

1 ^ I The Krawtchouk ensemble occurs in connection with zig-zag paths in random 

C^ ■ domino tilings of the Aztec diamond, and also in a certain simplified directed 

H ■ first-passage percolation model. We use the Charlier ensemble to investigate 

I the asymptotics of weakly increasing subsequences in random words and to 

cn ' prove a conjecture of Tracy and Widom. As a limit of the Meixner ensemble 

J^ ' or the Charlier ensemble we obtain the Plancherel measure on partitions, and 

fSj . using this we prove a conjecture of Baik, Deift and Johansson that under the 

^~~' ' Plancherel measure, the distribution of the lengths of the first k rows in the 

Q ' partition, appropriately scaled, converges to the asymptotic joint distribution 

On ■ for the k largest eigenvalues of a random matrix from the Gaussian Unitary 

^^ , Ensemble. In this problem a certain discrete kernel, which we call the discrete 

1 -Q I Bessel kernel, plays an important role. 

C^ ■ 

H ■ !• Introduction and results 

>'■ 

During the last years there has been a lot of activity around the problem of 

5^ ■ the distribution of the length of a longest increasing subsequence of a random 



X 



C^ 



permutation, its generalizations and their connection with random matrices, 
see for example [Ge], [Ra], [BDJl], [Jo3], [Ok], [BR2], [Bi], and also [AD] for 
connections with patience and the history of the problem. Let vr be a random 
permutation from the symmetric group Sn with uniform distribution Pperm,Af 
and let -L(vr) denote the length of a longest increasing subsequence in vr. It is 
proved by Baik, Deift and Johansson in [BDJl] that 

(1.1) lim Pperm,iv[^(7r) < 2Vn + tN^^^] = F{t), 
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where F{t) is the Tracy-Widom distribution, (1.5) for the appropriately scaled 
largest eigenvalue of a random M x M matrix from the Gaussian Unitary 
Ensemble (GUE) in the limit M ^ oo, see [TWl]. The probability density 
function on M for the M eigenvalues xi, . . . , xm of an M x M GUE matrix 
is 

1 ^ 

(1-2) 4>G\JE,Mix) = -TT- IT (2^i - ^of IT e'""' ' 

l<i<j<M j=l 

where Zm = (27r)^'^2^*^ '^ ni=i(iO "'^- This probability density can be ana- 
lyzed using the Hermite polynomials, which are orthogonal with respect to the 
weight exp(— x^) occurring in (1.2). Using standard techniques from random 
matrix theory, see [Me] or [TW2], we can write 



(1.3) Pgue.m 
where 



/ i 

max Xfc < V 2M + 



i<k<M ~ V2M1/6 



det(/ - /Cm) 



\L'2{t,oo)' 



ICM(tv) = ^^ , Km (V2M+ ^^ , ,V2M+ ^^ , | . 
Here Km is the Hermite kernel, 

T^ , ^ KM~lhM{x)hM-l{y) - hM-l{x)hM{y) -(a;2+«2)/2 

KM{x,y) = e ^ ^y -'/ 

KM x-y 

with hm{x) = Kmx"^ + . . . , /jg hn{x)hm{x) exp(— x^)(ix = (5„m; the normalized 
Hermite polynomials. It follows from standard asymptotic results for Hermite 
polynomials that 

. X , ,^ . X .. s Ai(nAi'(r/)-Ai'(OAi(r/) 
1.4 hm /Cm e,^ = ^ C,^ = —^ —, ^^^iZ, 

the Airy kernel, and also that the Fredholm determinant in the right-hand side 
of (1.3) converges to 

(1.5) F(t) = det(/ - A)L.(,^^) = j; ^ / det[A(e„e,)]t=i^'?' 

the Tracy-Widom distribution. 

The problem of the length of the longest increasing subsequence in a ran- 
dom permutation is closely related to the so called Plancherel measure on 
partitions, which occurs as a natural probability measure on the set of all 
equivalence classes of irreducible representations of the symmetric group. Let 
A = (Ai, A2, . . . , A£,0, 0, . . . ), Ai > A2 > • • • > Xi > 1, J2j ^j = A^, be a parti- 
tion of N, which can be represented in the usual way by a Young diagram with 
i rows and Xj boxes in the j row, see e.g. [Sa], [Fu]. Let / be the number 
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of standard Young tableaux of shape A. The Plancherel measure assigns to A 
the probabihty 

(1-6) Ppian,iv[{A}] = ^. 

The probabihty measure (1.6) is the push-forward of the uniform distribution 
on Sn by the Robinson-Schensted-Knuth (RSK)-correspondence, see e.g. [Sa] 
or [Fu], which maps a permutation tt to a pair of standard Young tableaux of 
the same shape A, and the length Ai of the first row is equal to L{tt). Thus, 
the length of the first row behaves in the limit as A^ ^ oo, as the largest 
eigenvalue of a GUE matrix. It was proved in [BDJ2] that the distribution of 
the rescaled length of the second row, Ppian,Af[A2 < 2VN + tN^'^], converges 
to the Tracy- Widom distribution for the second largest eigenvalue of a GUE 
matrix, [TW2], and it was conjectured that the analogous result holds for the 
k row. This conjecture will be proved in the present paper. It has recently 
been independently proved by Borodin, Okounkov and Olshanski, [BOO], see 
below. The conjecture also follows from the result by Okounkov in [Ok]. His 
proof uses interesting geometric/combinatorial methods. There are many ear- 
lier indications of connections between the Plancherel measure and random 
matrices for instance in the work of Regev, [Re], and Kerov, [Kel], [Ke2]. 

Another measure on partitions, coming from pairs of semi-standard 
tableaux, arises in [Jo3], where a certain random growth model is investigated. 
This measure relates to a discrete Coulomb gas on N of the form 

1 ^ 

(1-'^) Y~ n (.h,-h,fllw{h,), /iGN^, 

^ i<i<j<M j=l 

where Zm is a normalization constant. The weight w{x) = {^\^~ )q^, is the 
weight function on N for the Meixner polynomials, m„''^(x), see [NSU]. This 
measure on N^ can be analyzed using the Meixner kernel 

(1.8) K^lMi^,y) 

_ -q mM{x)mM~i{y) - mM-i{x)mM{y) ^ , . , xxi/2 

~ (\ -aVP x-v (w(x)w[y)) , 

with dn = n\{n + K — 1)!(1 — q)^ q^"'[{K — 1)!]^"*^, in much the same way as 
(1.2) is analyzed using the Hermite kernel. The Meixner kernel occurs in con- 
nection with probability measures on partitions also in the work of Borodin and 
Olshanski, [BOl]. The connection between certain measures on partitions and 
discrete Coulomb gases with their associated orthogonal polynomials is central 
in the present paper, and give them a very interesting statistical mechanical 
interpretation very similar to Dyson's Coulomb gas picture of the eigenvalues 
of random matrices. The difference is that in (1.7) we have a Coulomb gas on 
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the integer lattice instead of on the real line. Other statistical mechanical as- 
pects of measures on partitions have been investigated by Vershik, see [Ve] and 
references therein. We will refer to (1.7) as a discrete orthogonal polynomial 
ensemble. We will also be concerned with the cases w{x) = a^e~°'/x\, x G N, 
the Charlier ensemble, w{x) = {^)p^q'^~^, x G {0, ...,A^}, the Krawtchouk 
ensemble and w{x) given by (5.19), the Hahn ensemble. 
Consider the Poissonized Plancherel measure, 

(1.9) PglanKA}] = e-" Y. Ppian,iv[{A}]^, 

N=0 

on the set of all partitions, IPpian,Af[{A}] = if ^ • Xj / N. We will prove that 
this measure is a limit as g ^ of the Meixner ensemble. The Meixner kernel 
(1.8) converges in this limit, {q = a/M"^, K = 1, M ^ oo), to the discrete 
Bessel kernel 

(I 10) B"(x v) = /^ ■^^(^V^)'^;/+i(2\/q) - Jx+i{'2V^)Jyi2y/a) 

x-y 

This result can be used to give a new proof of (1.1), and also to verify the 
k row conjecture of [BDJ2], as well as to obtain asymptotic results in the 
"bulk" of the Young diagram. These results have recently been independently 
obtained by Borodin, Okounkov and Olshanski, [BOO], as a limiting case of 
the results in [BOl]. See the paper [B02] for a discussion of the connections 
between [BOO] and the present paper. 

The results for the Poissonized Plancherel measure can also be obtained 
as a limit of the Charlier ensemble. This ensemble arises in the problem of the 
distribution of the length of a longest weakly increasing subsequence in a ran- 
dom word which will be studied below. The random word problem has recently 
been investigated by Tracy and Widom, [TW3], using Toeplitz determinants 
and Painleve equations, see also [AD]. 

Before stating our results precisely we must introduce some notation. Let 

p = {A G N^+ ; Ai > As > . . . and ^ \j < cx)} 

j 

denote the set of all partitions, and p(^) = {A G P ; ^j Xj = N}, N > 0, 
the set of all partitions of N. Set ^(A) = maxjfc; Xk > 0}, the length of A. 
We will consider functions on V of the following form. Let / : Z ^ C be a 
bounded function which satisfies f{n) = 1 if n < 0. For a given L > we 
define g : V ^ Chj 



(1.11) giX)=llfiK + L 



i=l 
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We say that g is generated by /. Let Ql denote the set of ah functions g 
obtained in this way and write c{g) = ||/||oo- Let Vm = {A G "P ; i{X) < M} 
and pjf ^ = Pm n P(^). We also define, for M > 1, iV > 0, 

^M = {A G N^ ; Ai > As > • • • > A^}, 

M 

oif = {AGf]M; J^A, = iV}. 

Note that there is a natural bijection between Vm and Qm (and V\^ and 
n^^^). If M > L, g e Gl and X e Vm, then 

M 

(1.12) g{X)=l[f{X, + L-i), 

1=1 

since f{n) = 1 if n < 0, and we take (1-12) as our definition of g on ^m- 
For 771 > 1 and X £ V we define 

K^(A) = n (Ai-Aj+j-i), 

l<i<j<m 

and 



Wm{x) = n 



L (Aj + m- i)! 

According to a formula of Frobenius, see e.g. [Sa] or [Fu], the quantity /^ 
above can be expressed as 

(1-13) /^ = iV!^^(A)(A)H^^(A)(A). 

Let q £ (0, 1) and N > M. We define the Meixner ensemble on Qm by 

(L14) 



M-l ,,, ,,x, M 

,. ^\^^^ 

■Me,M,Nl 



[/ATI n ..^MivfT (iV-M)! - /A, + Ar-A 

[{A}] = (1 - g) 11 ^.!(^_m + ^-),^m(A) 11 1^;^^ + ^ _ ,J' 



Note that if we make the change of variables hi = Xi + M — i this gives us the 
discrete Coulomb gas (1.7) with the Meixner weight w{x) = {^\^~ )q^, where 
K = N — M + 1. For more about the Meixner ensemble and its probabilistic 
interpretations see [Jo3]. We can now state our first theorem. 

Theorem 1.1. For any g £ Gl, L >0, and a > we have that 
(1-15) E^ianb] = lim <S,7vb]- 

Thus the Poissonized Plancherel measure can be obtained as a limit of the 
Meixner ensemble. The theorem will be proved in Section 2. 
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Next, we define the Charlier ensemble on 0,m, which can be obtained as 
a limit of the Meixner ensemble, see (3.1). Given a > we define 



/ M-l \ M 

(1.16) PSh,M[{A}] =117? yMWWM{X)ll 



M-l , \ M 

■ ^ a \^i 
e 



3=1 -^ / 1=1 



a\^^ _^IM 



M 



on f^M- Again, the change of variables hi = Aj+M — i gives a discrete Coulomb 
gas, (1.7). The Poissonized Plancherel measure can also be obtained as a limit 
of the Charlier ensemble. 

Theorem 1.2. For any g £ Gl, L>0, and a > 0, 

(1-17) nUa] = Jim E:°ch,Mb]- 

M^oo 

The Charlier ensemble has a probabilistic interpretation in terms of ran- 
dom words, see Proposition 1.5. Since the Meixner and Charlier ensembles 
both correspond to discrete orthogonal polynomial ensembles they can be an- 
alyzed in a way similar to that in which the Hermite ensemble (CUE) is an- 
alyzed. This makes it possible to prove the following theorem, compare with 
[BOO]. 

Theorem 1.3. Let g G Ql, L > 0, be generated by f, see (1.11), and 
write (p = f — 1. Then, 

oo ^ k 

(1.18) E-pi,„b] = E M E n '^(^^■) det[S"(/i. - L, h, - L)]l^^„ 

k=0 ' heN*^ j=i 

where B"' is the discrete Bessel kernel, (1.10). Note that the right-hand side 
is the Fredholm determinant of the operator on £^(N) with kernel B°^{x — L, 
y-L)(j){y). 

The theorem will be proved in Section 3. 

As an example we can take (p{t) = —X(n,oo)i^) ^^"^ L = 0. This gives 

nUXi<n] = detiI-B-)l,^^^^^_^^^^^^^y 

By Gessel's formula the left-hand side is also a certain Toeplitz determinant, 
see e.g. [BDJl], and hence we get an interesting identity between a Toeplitz 
determinant and a certain Fredholm determinant on a discrete space. This 
formula has recently been generalized by Borodin and Okounkov, [BoOk]. 

By letting a go to infinity we can use (1.18) combined with de-Poisson- 
ization techniques to prove asymptotic properties of the Plancherel measure. In 
particular the next theorem generalizes the results of [BDJl] and [BDJ2]. Note, 
however, that we do not prove convergence of moments of the appropriately 
rescaled random variables. In Section 3 we will prove 
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Theorem 1.4. Letx^^' denote the j largest eigenvalue among the eigen- 
values xi, . . . , xm of a random M x M matrix from GUE with measure (1.2). 
There is a distribution function F{ti, . . . ,tk) on M^, see (3.48), such that 



(1.19) lim Pgue.m 

Af— +00 



x^^') < V2M+ ^ ^ , ,?• = 1, 



F(ti,...,tfe), 



for (ti,. . . ,tfc) G M'^, and 



(1.20) lim Ppian,7v[Aj <2^/iV + tjiVi/6,j = l,...,A;]=F(ti,...,tfc), 

for{ti,... ,tk)£R''. 



We turn now to the random word problem. By a word of length iV on M 
letters, M, A^ > 1, we mean a map w : {1, . . . , N} -^ {1, . . . , M}. Let Wm,n 
denote the set of all such words, and let Pw,M,7v['] be the uniform probability 
distribution on Wm,n where all M-^ words have the same probability. A 
weakly increasing subsequence of tf is a subsequence wi^ii), . . . ,w{im) such 
that ii < ■ ■ ■ < i„L and w{ii) < • • • < w{im)- Let L{w) be the length of a 
longest weakly increasing subsequence in w. The RSK-correspondence defines 
a bijection from Wm,n to the set of all pairs of Young tableaux (P, Q) of the 
same shape A G pv")^ where P is semistandard with elements in {1, . . . , M} 
and Q is standard with elements in {1,...,A^}. Under this correspondence 
L{w) = Ai, the length of the first row. Note that we must have £(A) < M, 
so £ G Vj^ which we can identify with Qj^ . In this way we get a map 
S:WM,N^dM\ 

Proposition 1.5. The push-forward of the uniform distribution on Wm,n 
by the map S : Wm,n -^ ^m ^^ 

(1.21) Pw,M,iv[5"'(A)] = Pch,M,N[{A}] = ^^( n -^)VMiXfWMiX) 

i=i ■^' 

on ^]^ ■ The Poissonization of this measure is the Charlier ensemble (1.16). 
Consequently., 

(L22) ¥^^mAL{w) <t\= Pch,M,N[Ai < t], 

and for the Poissonized word problem, 

oo N 

(1.23) P^,m[^(^) < i] = E ^"'^^^mMLH <t]= PSh,M[Ai < t]. 

N=0 
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Proof. See Section 4. 



D 



The probability (1.23) can also be expressed as a Toeplitz determinant 
using Gessel's formula, [Ge], see also [TW3] and [BRl]. The formula (1-21) can 
be used to prove a conjecture by Tracy and Widom, [TW3]. This conjecture 
says that the Poissonized measure on il,M induced by the uniform distribution 
on words converges, after appropriate rescaling, to the M x M CUE measure 
(1.2). In Section 4 we will prove 



Theorem 1.6. Let g be a continuous function onM . Then 

n OA^ V w \ f ^i-N/M Xm-N/M \ 

(1.24) lim Ech MN gi — ,. . . , , 

^ ^ N^^ ' • r V \/2iV7M y^m/M J_ 

= MIVttM / g{x)(t>GVE,Mi^)dxi ■ ■ ■ dxM-i, 

JAm 
where Am = {x G M^ ; zi > • • • > xm and xi + ■ ■ ■ + xm = 0}. Furthermore 

(1.25) hm Ech,M 



= M\ I 

J{xGR*f ;zi>-->a;M} 



g{x)(f)GVE,Mi^)d X. 



The case when g only depends on Ai has been proved in [TW3] using very 
different methods. 

The formula (1.23) can be used to analyze the asymptotics of the random 
variable L{w) on Wm,n as both M and N go to infinity. 

Theorem 1.7. Let F{t) be the Tracy-Widoni distribution function (1.5). 
Then, for all t £ M, 



(1.26) lim P^ 



M 



L{w) < ^ + 2 V^ + ( 1 + ^^ '"\^IH 



a 
M 



M J 



Fit). 



Assume that M = M{N) ^ co as N ^ co in such a way that (log Ny/^/M{N) 
-^ 0. Then, for all t G M, 



(1.27) lim Pw,M,7V 



LHS^ + 2V^+(l + f)^'V. 



F{t). 



Proof. See Section 4. 



D 
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Note that when M ^ a, the leading order of the mean goes hke 2y/a 
and the standard deviation hke a^'^ just as for random permutations. When 
M <^ a, we expect from (1.3) and (1.25) that 

L{w) = Ai « a/M + A/2a/M(\/2M + t/V2M^^^) 
= a/M + 2^/^ + t^jM^I'^, 
which fits perfectly with (1.26). 

In Section 5 we will consider two problems in combinatorial probability 
that relate to the Krawtchouk ensemble, namely Seppalainen's simplified model 
of directed first-passage percolation and zig-zag paths in random domino tilings 
of the Aztec diamond introduced by Elkies, Kuperberg, Larsen and Propp. 
Since both problems require some definitions we will not state the results here. 
A third problem, random tilings of a hexagon by rhombi, which is related to 
the Hahn ensemble will also be discussed briefly. 



2. The Plancherel measure as a limit of the Meixner ensemble 

The setting is the same as in [Jo3]. Let M.f^ denote the set of all A'' x A'' 
matrices with elements in N. We define a probability measure, P^[-] on A^at 
by letting each element Ojj in A G A^at be geometrically distributed with 
parameter q G (0, 1), and requiring all elements to be independent. Then 

(2.1) ^l\A\ = {\-qf\^^^\ 

A G A^AT, where S(^) = X^j j=iOij- Let A^Af(fc) denote the set of all A in 
A^AT for which ^(A) = k. Note that by (2.1) all matrices in A4i\i{k) have 
the same probability. Furthermore we let A4]sf{k) be the set of all matrices 
A in A^Af(^) for which "^iCHj < 1 for each j and '^jdij < 1 for each i; 
Mn = UfcA^Ar(fc). By taking the appropriate submatrix of A G Mi\f{k) we 
get a permutation matrix and hence a unique permutation. This defines a 
map R : M.j^[k) -^ Sk, where Sk is the k^^ symmetric group. Note that if q is 
very small a typical element in A^at belongs to Ai^^k) for some k. This is the 
crucial observation for what follows. The RSK-correspondence defines a map 
K : Mj\j{k) -^ 'py^i ^ and also a map S : Sk ^ 'pw , The number of elements 
in Sk that are mapped to the same A equals (/'^)^. It is not difficult to see 
that if A G Mnik) then K{A) = S{R{A)). Let g G Ql- It is proved in [Jo3] 
that 

(2.2) E^[<7(K(A))]=E^^_^^^[<7(A)]. 
With these preparations we are ready for the 
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Proof of Theorem 1.1. By (2.2) we see that in order to prove (1.15) it 
suffices to show that 



(2.3) lim E''J''\g{K{A))]=E'^pU9]. 



Note that F%[{A}\i:{A) = k] = 1/#7W^(A;), where #MN{k) = (^'~^+™), 
and F%[^{A) = k] = #7Wiv(A;)(l - q)^\^, by (2.1). Thus 

(2.4) 

oo 

E%[g{K{A))xMjA)] = J2K[9iK{A))xMjA)\j:{A) = k]F%[J:{A) = k] 

k=0 

oo 

= E E 5(^(^))(i - ?)"^V 

fc=OyleMiv(fc) 

oo 

= {l-qf'J2^' E 9{m{AeMN{k);K{A)=X}. 

k=0 AePC') 



The number of matrices in M]\f{k) which are mapped to the same permutation 
by i? is ( ^ ) , since there are ( ^ ) ways of choosing the rows and ( ^ ) ways of 
choosing the columns that select the submatrix. Since K = S o R we obtain 



#{A G M^ik) ; K{A) = A} = [^\ {f^f . 



Together with (2.4) this yields 



oo 



i%m{A))x^jA)] = {i-qf EfrTF^fcv^ E aW^- 

fc=o ■ ^ ^' xev(k) 



oo 



-fc AT\2 



= (1 - qf Y. l! (iV^l)!2^Pl-'^[^] 

CO J. ^ AT| \ 2 

= (l-a/A.r'|;^(^,j^-^)Ep,.„,bl, 

if we pick q = a/N"^. Since N\{N {N — k)\)^^ < 1 and converges to 1 as 
N ^ oo for each fixed k and furthermore Epian.fcis'] ^ c{g)^'^^^ ' ' , it follows 
from the dominated convergence theorem that 

(2.5) lim E^/^'b(K(A))x^,(^)] = EVnb]. 
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To deduce (2.3) from (2.5) we have to show that if M%- = Mn \ Mn-, then 



(2.6) 



,a/iV 



«'^^ [g{K{A))xM {A)]=Q. 



By the Cauchy-Schwarz' inequaUty, the left-hand side of (2.6) is 

(2.7) < ^''f\9{K{A)nl^FT\M),]"\ 

If A = K{A), then l{\) < S(^) and from the definition (1.11) of g it follows 
that 

\g{K{A))\ < c{g)'^^^(L,m) < c{g)^+^^^l 



Thus, 



Since, 



E-J^'^giKiA))'] < cigrY.<9r^T\n^) = k]- 



fc=o 






Af2 



iV^-l + A;\ / ayfaY 



[T.{A) = k] : 

as A^ ^ oo, it is not hard to show that 

(2.8) ¥!'f\g{K{A)f]<C{a,g), 

for all A^ > 1, where C(a,g) depends only on a and c{g). 
Next, we note that 



k\ 



M%c 



i=l y j J j=l I i J 



and hence 



tP'I^'^Xka* 1 ^ o^rTo"/^^ 



I^AT [A^^] < 2Arp' 



AT 



^ Ojj > 2 



Since, P?^[Ej %' > 2] = 1 - (1 - g)^ - iV(l - g)^-\ we obtain 

Ca2 



F^ [-A^^] < 



A^ 



Together with (2.7) and (2.8) this implies (2.6) and we are done. 



D 



It is also possible to give a more direct proof based on the explicit formulas 
similarly to what will be done with the Charlier ensemble in the next section. 
Above we have emphasized the probabilistic and geometric picture. 
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3. The Plancherel measure as a limit of the Charlier ensemble 

3.1. The limit of the Charlier ensemble. The Charlier ensemble is defined 
by (1.16). It can be obtained as a limit of the Meixner ensemble (1.14) by 
taking q = a/MN and letting N ^ oo with M fixed. In this limit 

(3 1) (1 a)^^W ^^"'^^' ftP- + ^"^^/' 



n 71^^(^)11 



r^V"e-°/M 



VM/ 



so we obtain (1.16). In light of Theorem 1.1 we see that it is reasonable to 
expect that the Poissonized Plancherel measure should be the limit of the 
Charlier ensemble as M ^ oo. The interpretation of the Charlier ensemble 
in connection with random words, Proposition 1.5, also supports this since a 
random word in the limit M ^ oo is like a permutation (no letter is used 
twice), see also [TW3]. We will give an analytical proof of Theorem 1.2 that 
does not use the RSK-correspondence. We start with the following simple but 
important lemma. 

Lemma 3.1. If M > £{X), then 

(3.2) Vm(A)W^m(A) = 1^^(A)(A)VF,(A)(A). 

Proof. We may assume that M > i{X). Note that, by definition, Aj = 
if i > ^(A). Hence, 

iW M 

^M(A)=y,(A)(A)n n (A.+i-i) n ^^-^^ 

i=i j=e(x)+i e{\)<i<j<M 

T. /xxTT (Aj + M-i)! TT 

i=l ^ * ^ ' ' e{X)<i<j<M 

Thus in order to prove (3.2) we must show that 

M 

n (j'-^)= n (A.+M-0!, 

e{x)<i<j<M i=e(x)+i 

but this is immediate since Aj = if i > ^(A). D 

Proof of Theorem 1.2. It follows from the definition (1.12) of g{X) that 

(3.3) |<7(A)| < c(g)"^^^W^)'^) < c((7)^W+^. 
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Let Pch,M,N be defined by (1.21). Then, 

(3.4) 

/m-1 \ 
Ech,Mb] = E 5(A) n 7i VMiXfWMWH 



M 



Aen 



M 



i=i 



i=i 






E«""^ E 5(A)Fch,M,N[{A}]. 



Ar=0 



Nl 



AeQ 



(iV) 



Thus, by (3.3) and the fact that £(A) < A^ if A G n^^\ 



(3.5) 



j; 5(A)Pch,M,N[{A}] 



Aen 



(AT) 



<c(<7) 



L+Af 



)W 



since Pch,M,N is a probabihty measure on ^^^ ■ Given e > we can choose K 
so large that 



(3.6) 
Consequently. 

(3.7) 



°^ ^Af 



L+A 



N=K+1 



m 



<e. 



K N 

Ech,Mb] - E «""^ E 5(A)Pch,M,N[{A}] 



A=0 



Aeo 



(JV) 



<e. 



If M > K > Af > ^(A), A G 0S^J^\ we can identify ^-^^ with p(^) and 
(3.2) to write 

(3.8) Y. 5(A)Pch,M,N[{A}] 



use 



Aen 



(iV) 



M ., ... M-1 

E .(A)iv!^.(A)(A)^T^.(A)(A)^n ^^^ ^r n^ 



AeQ 



(iV) 



i=i 



i=i 



^(A) 



(A,- + M-i)! 



E 5(A)^Plan,A[{A}]n^X7M 



■JJ 



where the last equality is a straightforward computation using the fact that 
Aj- = Oifi >^(A). Now, 



iiM^^{M-j)\ 111 



j-iA ri_i-^(A) 



M 



M 
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which goes to 1 as M ^ co for a fixed A. Since the sum in (3.7) is, for a fixed 
K, a sum over finitely many A, we obtain 

K jv 

(3.9) Jim E^""^ E 5(A)Pch,M,N[{A}] 

= E«""^ E 5(A)Ppian,iv[{A}]. 

Using (3.6) and the fact that Ppian,7V is a probabihty measure on V^^\ we 
obtain 



(3.10) 



K AT 

^Planb] - E ^""^ E 5(A)Ppian,iv[{A}] 
Af=0 ' \^oW 



< £. 



The theorem now follows from (3.7), (3.9) and (3.10). D 

3.2. Coulomb gas interpretation of the Plancherel measure. As M ^ oo 
the number of particles in the Coulomb gas representation of the Charlier 
ensemble goes to infinity, so a Coulomb gas interpretation of the Plancherel 
measure is not immediate. We will now show that we can actually approximate 
'^Plan ^y ^ Coulomb gas with K particles, which gives a good approximation 
if K is chosen large enough (depending on a). 

Consider the Poissonization of the restriction of the Plancherel measure 
to V^^^ 



m9] = e--Y.m E aW 



m ^^ ^'^ ' m 



^=0 AepW 
for g G Ql. If M > L it follows from (1.12), (1.13) and Lemma 3.1 that 



M 

When M is large, we expect that -^^[5] and Epjjjjj[(7] should be close. 

Lemma 3.2. Assum,e that g G Ql and let d > be given. There is a 
numerical constant C such that if M > max(I/, aexp((i + 1)), then 

(3.11) |K"p„,b] - FU9]\ < C{c{g)e-'')^ . 

Proof. Set 
Rn,m[9] = Epianb] - FS,[g] = Y. e""— Y. 5(^)^- 

N=0 ■ AeP(iV)\p(f) 
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U N < M, then Rn,m[9] = since then £(A) <^iXi = N < M,so p(^) = 
V^^K UN > M>L, then \g{X)\ < c{g)^ since Ai = if i > iV. Thus, 



oo jv 

..a 



\ni.n[9] - FUaW < E ^~"j^\Rnm[9]\ 



N=M+1 



i E ^-"^^w" 



N=M+1 

This last sum is estimated as fohows. By Stirhng's formula there is a nu- 
merical constant C such that exp{—a)a^/N\ < C exp{—af{N/a)), where 
f{x) = xlogx + 1 - X. If N/a > exp(<i + 1), then f{N/a) > dN/a, and 
so exp[—a)a /N\ < ex.p{—dN). The lemma is proved. D 

Recall from the introduction that Vm can be naturally identified with 
Qm- For K < M we define 

^M,K = {A G Qm ; Xk+1 = ■ ■ ■ = Am = 0}, 

and Q*j^ ^ = Qm \ ^m,k ■ If 1 < j < M — K we set 

^*m,kU) = i^ ^ ^*M,K ; Am+i-j > but Ai = 0, M + 1 - j < i < M}, 

so that ^*M K ~ ^j=i ^*M kU)- '^^^ next lemma asserts that ^(A) is not too 
large for typical A that we will consider. 

Lemma 3.3. Let g G Gl be generated by f. Assume that f satisfies 

(3.12) < /(x) < Co/(x - 1) 
for all X £ Tj and some constant Cq. Then 

(3.13) e'" Y. aWVMixfWMW'Ha'^ < (^'_ ■ . ^..^ ^Mia]- 

Proof. It will be most convenient to use the discrete Coulomb gas rep- 
resentation. Set Xj = Xm+1-j + J — 1, J = 1,...,M and let Am{x) = 
ni<i<7<M(^i ~ ^«) ^^ ^^^ Vandermonde determinant. Also, set A = {x £ 
N^ ; < xi < • • • < xm} and Aj = {x G ^ ; Xj < i for i < j and Xj > j}, 
j = 1,...,M. Note that A G ^*m kU) translates into x G Aj. If x G Aj, 
then Xi = i — 1 for i = 1, . . . ,j — 1 and we have the first hole in the particle 
configuration x at j — 1. Now, 
(3.14) 

Y. 9{X)VM{XfWM{Xflla'^^ = Y AM{xfll^f{x^ + K-M). 
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We want to show that, with high probability, the first hole must be fairly close 
to M. Define Tj : Aj ^ A by Tj{x) = (xi, . . . , Xj-i,Xj — 1, ... , xm — 1) = x' . 
Clearly, Tj : Aj -^ Tj{Aj) is a bijection. Write 

M ^. 

L%{x) = AM{xf n ^/(^^ + K-M). 
1=1 

For X £ Aj, 

2 M , ,,so / . / N \ 2 M 



^Mix')J l{ix,\r \AMix')J 11 Xf 



Since 



and 



n-<TT- = f^^^^^ 

11x2-11^2 I M! 

Am{x) _ T-r Xfc ^ T-r A: _ f M 

^m{x') ~ 11 xfc - (j - 1) - 11 k-{j-l)~ \j - 1 
if X G Aj, we obtain, using our assumption on /, 

Inserting this into (3.14) yields 

M 
e-" j; 5(A)yM(A)2|^Af (A)2 J] a^^ = e"" ^ L^(x) 

^(M-i + l)!2" J^^m(T,(x))<^^_^.^^^,,F^M, 

and the lemma is proved. D 

Lemma 3.4. Let g £ Ql be generated by f which satisfies (3.12). Assume 
that M > K > max(L, e\/2Coa). Then, 

(3.15) \F-[g]-F^[g]\<2(^j^^^ ^ F^b]- 

Proof. If A G r^M.i^ then ^(A) < K < M and hence by Lemma 3.1, (1.11) 
and the fact that ^m,k and Qk can be identified we obtain 

M 

e- Yl 9WVM{XfWM{Xflla^^ = F^[g]. 
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The left-hand side of (3.15) is 

M-K M 

^ E ^~" E 5(A)VM(A)2w^M(A)2n«'' 

by (3.13). This completes the proof of the lemma. D 

We can now demonstrate how the Plancherel measure can be approxi- 
mated by a Coulomb gas, (compare with the discussion in the Appendix in 
[BDJl]). 

Proposition 3.5. Assume that g G Qk is generated by f which satisfies 
(3.12). Let K = [r^, r > V2Q^. Then, 

(3.16) E-pUg] = (l + 0(r-^))^ Yl ^^W' 11 r^ n/(^^)' 
where 

Proof. Write 

(3.17) E^Ua] - Epianb] - F^[g] + E^p^^Jl] F^b] F^W F^[l] ' 
By Lemma 3.2 

(3.18) hm FSi[g]=EWn[9]- 

M^oo 

By Lemma 3.4 and the choice of K and r 

(3.1.) i^ - + «<■-''■)■ 

for any M > K, and similarly with g replaced by 1. Using (3.18) and (3.19) 
in (3.17) and letting M ^ oo we obtain 



Ef.lanb] = (l + 0(r--")) 



X^^^xb] 



which is exactly (3.16). The proposition is proved. D 
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Thus we have an approximate Coulomb gas picture of the (shifted) rows 
of A under the Plancherel measure analogous to Dyson's Coulomb gas picture 
for the eigenvalues of a random matrix. 

Remark 3.6. The confining potential for the discrete Coulomb gas in 
(3.16) is 



with limit 



V^[h,] = -^log{a'^/{W.f) 



lim V^[Kx] = 2[xlogx+ (logr — l)x] = V{x). 



We can now use general techniques for Coulomb gases, see e.g. [Jol], [Jo3], 
to deduce asymptotic distribution properties. The potential V has the (con- 
strained) equilibrium measure; compare with Section 2 in [Jo3], u{t)dt, where 

1, if < t < 1 - 2/r 

i - ;^ arcsin(§(t - 1)), if 1 - 2/r < t < 1 + 2/r 

0, if t > 1 + 2/r. 



u{t) 



Pick f{t) = exp{(f){t/[ry/a])) with (/> : M ^ M continuous, bounded together 
with its derivative and (f){t) = if t < 0. Then, 



giX) = [[expi^i rr7=\ ))• 



r^/a] 

If we pick r sufficiently large (depending on 0) we can use (3.19) and (3.20) to 
show that 

I rl+2/r 

(3.20) lim ^^.logE'^i,J<7(A)]= / ,^(t)n(t)dt. 

" - [r^\ Jo 



+ 00 



From the limit (3.20) it is possible to deduce Vershik and Kerov's 0-law for 
the asymptotic shape of the Young diagram, [VK], see also [AD], where an 
outline of the argument using the hook-integral is given. (The r-dependence 
in the formulas above goes away after appropriate rescaling.) From what has 
been said above we see that the i7-law is directly related to an equilibrium 
measure for a discrete Coulomb gas. Using the general results in [Jo3] we can 
use (3.16) to show upper- and lower-tail large deviation formulas for Ai (= the 
length of the longest increasing subsequence in a random permutation) under 
the Poissonized Plancherel measure. These formulas have been proved in [Sel], 
[Jo2] and [DZ] by other methods and we will not give the details of the new 
proof. 

3.3. Proof of Theorem 1.3. We will now use Theorem 1.2 to prove The- 
orem 1.3, but before we can do this we need certain asymptotic results for 
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Charlier polynomials and Bessel functions. Let 

Wa(x) = e~" — , X G N, a > 0. 
x! 

The normalized Charlier polynomials, Cn{x;a), n > are orthonormal on N 
with respect to this weight. The relevant value of the parameter a for us will 
be a = a/M, and we define the Charlier kernel 

(3.21) 

T^a . . _ /- cm(x; ^)cM-i{y; ^) - cm^iJx; fj)cM{y; fj) 

-"■Ch,Ml^>yj — V" 

X y 

for X 7^ y and 
(3.22) i^^h,M(^,^) 

= V^W^/Mix) 

The polynomials Cn{x;a/M), n > 0, have the generating function 

oo ^ ^ n/2 -, 

It follows from this formula that we have the following integral representations. 
If < r < ^/M, then 



a \ / a \ / a \ , / a 



H, \ ■■'; jj) «'-i l,^; m) - «'-■ \'=- m) =m (^i; jg 



,3.23) c„(.; ii) = ,/irj- /' e-^'"" fl + ^V ' .. 



— 7r 



and if ^/a/M < r, then 
(3.24) 



c fx- ^) - . /^^ r e- v^-" Tl + ^^) ' ^ iP 



for any x G M, where the powers are defined using the principal branch of the 
logarithm. 

We want to write the Charlier kernel in a form that will be convenient for 
later asymptotic analysis. Define, for a given r > 0, x G Z, 

AUx) = V^j^w^/m{x) (I + —J e-2^/^, 
D^^^(x a)-— r aire'')e^^^-''''^ f ^ + Mre^' V 1 
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Fm'\x^9) 



J y/a/M 

-tcx.r 



y/a- Ms 



^ + M 



-m'^^ 



if r > yJajM, and if r < ^/a/M, then F^^{x,g) = 0. Then, some computation 

shows that when x is an integer (the case we are interested in), 

(3.25) 

DM''{x,gi)D'^/'{y,92) - D'^''{x,g2)D'^[{y,gi) 



^Ch,M\^iV) 



A%{x)A%{y) 



x-y 



when X 7^ y, and 

(3.26) 

^Ch,M(a;>a;) = Alj{x)[D'l^\x,g2)D'^^ {x - 1,53) - D'^''{x,gi)D'^''{x - 1,54)] 

+ Alj{x)[F^''{x,gi)D'^''{x,g2)-F^''{x,g2)D'^''{x,gi)], 
where gi{z) = 1, 52(2) = ^ - 1, 

and 54(2) = g2{z)g^{z). Note that all the giS are bounded on |z| = r. 
The discrete Bessel kernel is defined by (1.10) for x ^ y and 

(3.27) B"(x,x) = V^[L,(2V^)J,+i(2V^) - J,(2V^)L,+i(2V^)] 

for x = y, where Lx{t) = -^Jx{t)- The Bessel function has the integral repre- 
sentation 

(3.28) 



^V^(^e-^''-re^'^)+ixe^x^Q 



smvrx r ^^(_i/s+,)^d£ 



71" JO 

for X G M, r > 0. Differentiation shows that for integer x, 

1 r'^ 



(3.29) 



Lx{2V^) 



2tt 



log{re'^)e'^^-^ '"-re^n+i^dr'^de 



-1 






/O s 

The next lemma shows that the discrete Bessel kernel is the M ^ 00 limit 
of the Charlier kernel and establishes some technical estimates. (We will only 
consider the case when x, y are integers but this restriction can be removed.) 



Lemma 3.7. For any x,y G 
(i) 
(3.30) 

(ii) 



hm K"ci,^j^{M + x,M + y) = B''{x,y). 



(3.31) 



B"(x,y) = ^ J,+fe(2V^)J^+fe(2V^). 



k=l 
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Furthermore, there is a constant C = C{a,L), such that 
(iii) 



(3.32) 



J2 K'^c^^m{M + x,M + x)<C 



if M is large enough, and 
(iv) 

(3.33) 



Y^ 5°(x,x) < C. 



{In (3.33) we can take C{a,L) = a/\/2 + L.) 

Proof. We have to show that (3.25) and (3.26) converge to (1.10) and 
(3.27) respectively. Using Stirhng's formula we see that A'^{M + x) ^ ^/a as 
M ^ oo. The result then follows from the integral formulas above, the fact 
that 






^Va{l/z~z)^y 



and gsiz) -^ zlogz as M ^ oo. This establishes (3.30). The identity (3.31) 
follows from the recursion relation Jx+i{t) = 2xJx{t)/t—Jx~i{t), which implies 

B^ix,y) = Jx+iit)Jy+iit) + B''ix + l,y + l), 

and (3.31) follows by using the decay properties of the Bessel function; see 
Lemma 3.9 below. 

The estimate (3.32) is proved using the formula (3.26). Stirling's formula 
can be used to show that A'^{x + M) < 2^/a for all x > 0. We have 



Va + Afz 



Va + M 



M 



\M 



\fa^Mz 



^/a + M 



< 1 + 



M 



Z\ + 



M\z\j y ' M 

<exp((l-(5)-iV^)(l-5/2)^ 



if |z| = r = 1 — 6 and M > 2y/a/6. This estimate can be used in the integral 
formulas for D.'r and F,) and we obtain 



'M 



M 



Thus, 



D^^/^M + x; g,)\, \Fli^'\M + x;^.)! < Ce^^i-f. 



J2 ^Ch,M(M + X, M + x) < CV^e'^ J2 (j)"- 



The estimate (iv) can be proved in a similar way but we can also proceed 
as follows. Using the generating function for the Bessel functions Jn(t), n G Z, 
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one can show that Yln=i n^Jnitf = t^/^, see [Wa, 2.72(3)], and Yln=i M^f = 
i(l - Mtf) < 1/2, so by (ii) and the fact that B{x,x) < 1, 

oo oo 

Y^ S"(x, x)<L + Y^ nJn{2^/af < a/V2 + L, 

x=—L n=l 

where we have used the Cauchy-Schwarz inequality. The lemma is proved. D 
We are now ready for the 
Proof of Theorem 1.3. We have 

M-l M A^ / \ A, 

i=i ■^' xeUM *=i *=i 

If we make the change of variables hi = Xi + M — i, this can be written 
^ M M 



'Ch,ML 



M /ig^M i=i i=i 



Now, using a standard computation from random matrix theory, see [Me], 
[TW2], we can write this as 

M-l k 

(3.34) E"ch,Mb] = E E X{4>{h,)dei{Kl^^M{K+M-L,h,+M-L)%^^ 

fc=0 ftgNfe «=1 

since 0(t) = if t < 0. The Charlier kernel is positive definite, so we have the 
estimate 

k 

I det(i^ch,M(2^J5^i))ij=i| ^ 11 -^ch,M(2^i'^i)• 
Thus, by Lemma 3.7(iii), 

k 

Y, n '^(^*) det(i^^h,M(^* + M - L, /i,- + M - L))f^^.^i 



Y i^^,M(M + x,M + x) <(C||</>||c 



The analogous estimate for the Bessel kernel follows from Lemma 3.7(iv). 
These estimates and Lemma 3.7(i) allow us to take the M ^ oo limit in 
(3.34). By Theorem 1.2 this gives (1.18). The theorem is proved. D 

Note that we could just as well use Theorem 1.1 and the Meixner ensemble 
to prove Theorem 1.3. The proof would be the same and we just have to prove 
(3.30) and (3.32) for the Meixner ensemble instead. 
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3.4. Asymptotics of the Plancherel measure. Theorem 1.3 can be used 
to analyze the asymptotic properties of the Plancherel measure in different 
regions. One can distinguish three cases corresponding to three different scaling 
limits of the Bessel kernel. First we have the edge scaling limit, 

(3.35) lim q1/6b"(2V^ + ^a^/^ 2^^ + 7]a^/^) = A{C, r?), 

a— »oo 

where A is the Airy kernel defined in (1.4). This is the case that is considered 
in Theorem 1.4. Secondly we have the bulk scaling limit, 

(3.36) lim B'^ir^, r^a + u) = — ^ -, 

It G Z, —2 < r < 2, where R = arccos(r/2); the right-hand side is the discrete 
sine kernel. We will not discuss the local behavior in the bulk of the Young 
diagram; see [BOO]. Thirdly we have an intermediate region, 

(3.37) lim 7rai/^-^/2^"(2V^ - a^ + i.ia'l^-^l^ ,2^ - a^ + T,T]a'l^'^l^) 

sin7r(.5 — r\) 
7r(^ - ri) 

if 1/6 < (5 < 1/2, the ordinary sine kernel. Thus in this region the local 
behavior is the same as that in the bulk in a random Hermitian matrix. The 
limits (3.35) to (3.37) can be proved using the saddle-point method on the 
integral formula for the Bessel function. From the point of view of the Coulomb 
gas picture of the Young diagram, the cases one and three are similar to the 
random matrix case since at the edge a discrete Coulomb gas approximates a 
continuous Coulomb gas. Case two is different however, since in the bulk the 
discrete nature is manifest; the charges sit close to each other. 

Before turning to the proof of Theorem 1.4 we have to say something 
about de-Poissonization, the joint distribution of the first k rows (A; largest 
eigenvalues) and the asymptotics of the Bessel kernel. 

We have the following generalization of a lemma in [ Jo2] . 



Lemma 3.8. Let ^^ = N + 4^/N log N and u^ = N - AyJNlogN. Then 
there is a constant C such that, for < Xi < N, 

(3.38) Pl^L[Ai <xi,...,Xk<Xk\- — < Ppia„,iv[Ai < xi,. . . ,Xk < Xk] 

<rpU>^,<x^,...,Xk<xk] + -^. 

Proof. This is proved as Lemmas 2.4 and 2.5 in [Jo2]. Denote a permu- 
tation in Sn by vr' -^ and let Sn+i{J) denote the set of all vr' ^^^ such that 
7]-(^+i)(A^ + 1) = j. Each vr' +^) in Sn+i{J) is mapped to a permutation 
Fj{7:^ ^^') in Sj\f by replacing each vr' ^^^(i) > j by tt^ ^^'{i) — 1. The map 



282 KURT JOHANSSON 

Fj is a bijection from Sn+i{J) to Sn- Apply the Robinson- Schensted corre- 
spondence to Fj{Tr^'^^^') to obtain the P-tableau. Replace the entries i by 
i + 1 for i = j, . . . ,N and then insert j. This insertion can only increase the 
length of any row and we obtain the P-tableau for vr' ^^^ . Thus, 

A.(F,(vr(^+i))) < A.(vr(^+i)), 

for all rows. If we define g{TT^^') to be 1 if Ai(7r^) < Xj for i = 1, . . . , fc and 
otherwise, we see that 

<7(F,(^(^+i)))>5(vr(^+i)), 

and we can proceed exactly as in [Jo2] using the fact that the Plancherel 
measure on V^^' is the push- forward of the uniform distribution on Sn. D 

For X G M , n G N and a sequence I = (/i, . . . , I^) of intervals in M we 
let x{'^:^j x) denote the characteristic function for the set of all x G M^ such 
that exactly Uj of the Xj's belong to Ij, j = 1, ... ,k. A computation shows 
that for a single interval 

1 a"i ^ 

X{lj,nj,x) = —^^-n-YlC^ + ZjXiiixi)) 

nj.az- .^^ 

and hence 

(3.39) X(^, ^' ^) = ^. ! . . . nk\ 9^"^ ■ . . 9z"^ H H^^ + ^^^h (^»)) 

1 ' ' ' 3 1=1 jf=l 



=2fc=-l 



Note that if the intervals are pairwise disjoint, then J|-^]^(l + ZjXi\xi)) = 
1 + Yl,j=iZjXij{xi)'-, compare with [TW2]. Let P be a probability measure 
on M*^ and let ai > ••• > a^. Set Ij+i = (0^-1-1,0^], j = 1, . . . ,k — 1 and 
/i = (ai, 00). Let 



Lfc = {n G N^ ; ^ nj < r - 1, r = 1, . . . , A:}. 
i=i 

Define x^^> to be the j^^ largest of the Xj's. Then, 

(3.40) P[x(^) < ai, . . . , x(^) < Cfc] = ^ E[x(X, n, x)]. 

Hence, the problem of investigating the distribution function in (3.40) reduces 
to investigating expectations of the right-hand side of (3.39). 

In the proof of Theorem 1.4 we will need some asymptotic results for 
Bessel functions. 
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Lemma 3.9. Let Mq > be given. Then there exists a constant C 
C(Mo) such that if we write x = 2y/a + ^a^'^ , then 



imin^y/Me|V2)|e| 



(3.41) \J^{2^/^)\<Ca~'^/^exp 

for ^ G [—Mq, oo). Furthermore 

(3.42) lim a^/V^(2Va) = Ai{C), 
uniformly for ^ G [—Mq, Mq]. 

This can be deduced from classical asymptotic results, [Wa] and it is also 
rather straightforward to proceed as in Section 5 of [Jo3] using the integral 
formula (3.28). 

We are now ready for the 

Proof of Theorem 1.4. We will prove (1.20). The proof of (1.19) is analo- 
gous using the Hermite kernel instead. From Lemma 3.8, the fact that a dis- 
tribution function is increasing in its arguments, that the distribution function 
F{ti, . . . ,tfc) is continuous and yJJxN — \fN ss 2y/\ogN, ,Ji^—\fN « 2^/log A^, 
we see that it suffices to prove that 
(3.43) 
lim PgiaJAi - 1 < 2^^ + tlal/^ ...\k-k<2^ + tka^'"^] = F(ti, . . . , t^), 

a— >oo 

for any fixed (ti, . . . , tfc) G IR^, ii > • • • > tfe- Set 

Ij+i = {2^ + tj+W^^,2^ + tja^l% j = l,...,k-l 

and /i = (2-^0-1- tia^'^, oo). By (3.39) and (3.40) it is enough to consider the 
expectations 

(3.44) E-pi,, 



oo k 



If we write 0a(s) = n7=i(l + ^jXli^)) — 1 it follows from Theorem 1.3, with 
L = 0, that the expectation (3.44) can be written as 

oo ^ k 

(3.45) FUz,t) = E^ E ll^a(.h,)det[B-{h,hj)]'l^^^,. 

k=o ' /ieN'=i=i 

Note that Fa{z,t) is an entire function of z. Set Jj+i = {tj+i,tj], j = 1,. . . , 
k — 1, Ji = (ti,oo) and write ip{s) = Y\j=iil + ZjXJj{s)) — 1. Define 



(3.46) 



oo ^ „ A; 

^(^'*) = Ery / n^(^^-)det[^(e„e,)]f,,=i. 

k=0 ■ ■^^'' 3=1 
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We want to show that 

(3.47) hm Fa{z,t) = F{z,t), 

a— »oo 

uniformly for z in a compact subset of C . Then also derivatives of Fa{z,t) 
converge to the corresponding derivatives of F{z,t). The limit (3.43) then 
follows with 

(3.48) Fit„...,t,)=Y: , ,g^n, Q.n^ n^^t) 

So it remains to prove (3.47). Note that 0q,(s) = if s < 2^/a+t|^a^'^ and 
that0a(s) = 'ip{a~^/^{s-2./a)). Given r G Mweset A(r) = {r,r+l,r+2, . . . }. 
Then, 

°" 1 ^ ' / /j, \ ^ ^^^ 1 



(3.49) FUz,t) = J2- J2 n^b7aVi-t[^"(^'-.aV6y' 

where B°'{^,r]) = a^^^B°'{2y/a + ^a^/^ ,2y/a + T]a^^^). We can now prove 
that (3.47) holds pointwise in z by the same argument as was used in the 
proof of the analogous statement in Section 3 of [Jo3] . That proof depends on 
the following properties of the kernel; compare with Lemma 3.1 in [Jo3] and 
Lemma 4.1 below. 
(i) For any Mq > there is a constant C = C{Mq) such that for all S, > —Mq 



Y, B^{2^ + ia^/^ + m, 2^ + ^a^/^ +m)<C. 

m=l 

(ii) For any e > 0, there is an L > such that 

oo 

Y^ B'^(2^f^ + La^l^ + m, 2^f^ + La^l^ + m) < e, 

m=l 

for all sufficiently large a. 
(iii) For any Mq > and any e > 



^a ( _^ ^\ _ ^ ^ " ^ 



<e 



for all integers m,n ^ [—Mqo^'^^Mqq^'^] provided a is sufficiently large. 

The estimate (i) is used to estimate the tail in the /c-summation in (3.49), 
(ii) is used to limit the /i-summation and (iii) is used to prove that the Riemann 
sums converge to integrals. 

If z belongs to a compact set K there is a constant C, independent of z, 
such that llV'lloo ^ C. Together with (i) this shows that the family {Fa{z,t)} 
is uniformly bounded for a > 0, z G -ftT and hence (3.47) holds uniformly by a 
normal family argument. 
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The properties (i) to (iii) above are straightforward to prove using the 
representation (3.31) and Lemma 3.9. To prove (i) and (ii) we use 

oo oo 

m=l n=l 

which can be estimated using (3.41) (we get a Riemann sum). Similarly, 
^"i~Tj6^ ^Te) can be written as a Riemann sum, using (3.31), which is con- 
trolled using (3.41) and (3.42). This Riemann sum can be compared with the 
corresponding Riemann sum for the following representation of the Airy kernel, 

[TWl], 



poo 

A{C,7])= A[{^ + t)Ai{7] + t)dt 
Jo 



'0 
and in this way we obtain (iii). 



4. Random words and the Charlier ensemble 

In this section we will prove our results on random words. 

Proof of Proposition 1.5. Let L{M,N, X) denote the number of pairs 
(P, Q) of tableaux of shape A G ^i]^^ with P semistandard with elements 
in {1, . . . , M} and Q standard with elements in {1, . . . , A^}. Then 

(4.1) ^wmMS^'W] = ]^^(^' ^' ^)- 

The number of possible P's is, by [Fu], 

l<i<j<M •' \i=l 

and the number of possible Q's is / given by (1.13). By (4.2), (4.3) and 
Lemma 3.1 we obtain 

(4.3) L(M, AT, A) = A! f [] -y J Vm{\?Wm{X)- 

Inserting the formula (4.3) into (4.1) yields the desired result (1.21). The 
formulas (1.22) and (1.23) are immediate consequences. The proposition is 
proved. D 

Next, we give the 



286 KURT JOHANSSON 

Proof of Theorem 1.6. We will prove (1.24); the proof of (1.25) is analo- 
gous. Both are straightforward asymptotic computations using Stirling's for- 
mula and we will indicate the main steps. Set 

Ao - N/M , ^ 

Note that X]7=i ^j — 0> since X^,=i Aj = N . Then, 

Af/M+M-i 



(\ , )l^ -M 2ttN ( N\ ' - x^-N/M+o(l) 



as N ^ oo, and hence 

-M/2 / ^^^ Af+M(M-l)/2 M_ 

e "'J. 
i=i 
Furthermore, 






and consequently 

(4.4) 

M-l M . . _(M-l)/2 

Pch,M,N[{A}] ~ V^(2vr)-^/22^^/2 n ^^M{xflle-^U^ 






= V^^M!0GUE,M(a:). 

From this we see that the left-hand side of (1.24) is approximately a Riemann 
sum for the right-hand side, which in the limit N ^ oo converges to the right- 
hand side. The factor M! in the last expression in (4.4) comes from the fact 
that in (4.4) the variables are ordered. This completes the proof. D 

For the proof of Theorem 1.7 we need asymptotic results for the Charlier 
kernel analogous to those for the Bessel kernel in the proof of Theorem 1.4. 

Lemma 4.1. Let u = M + a/M + 2y/a and cr = (1 + y/a/Mf/^a^/^ . 
(i) For any Mq > there is a constant C = C{Mq) such that, for all 
integers n > —M^a, 

oo 

(4.5) Yl ^Ch,M(M +n + m,[u]+n + m)<C. 

m=l 

(ii) For any e > there is an L > such that 

oo 

(4.6) Y. ^Ch,M(M + WL] + m, M + [aL] + m) < e 
if M, a are sufficiently large. 
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(iii) For any Mq > and any e > 0, 

(4.7) |^i^Sh,M(M + m, M + n) - ^ (^, ^) I < 5 

for all integers m,n £ [— MqO", Mqct] provided a and M are sufficiently large. 

Proof. The proof is based on the formulas (3.25) and (3.26) for the Char- 
her kernel. The proof is completely analogous to the proof of the corresponding 
result for the Meixner kernel in Lemma 3.2 in [Jo3, §5], so we will not give 
the details here. Asymptotic formulas for Charlier polynomials with fixed 
a = a/M have been obtained in [Go]. D 

Proof of Theorem 1.7. By (1.16) and (1.23) 

M M 

F'w,M[iH<s] = n- E AM{hfl[w^/M{h^), 

niaxhj<s+M~l 

where we have made the substitution hi = Xi + M — i. Using Lemma 4.1 
this can be analyzed exactly as the analogous problem involving the Meixner 
weight in Section 3 in [Jo3]. Lemma 3.1 in [Jo3] gives 



(4.8) 



W,M 



L{w)<^ + 2^+(l + ^\'^\^l^i 



M+^^"+^^ + 17; 



no, 



as a, M ^ cxD with -F(^) given by (1.5). This proves (1.26). Next, we observe 
that for fixed M, Pw,M,Af[^('U^) < s] is a decreasing function of N ^ which can 
be proved as the corresponding result for permutations in [Jo2]. Thus, with 
II N and vi\j as in Lemma 3.8, we have 

(4.9) P^>[L(«;) <s\--^< Pw,M,iv[i(u;) < s\ < P^,m[^(^) ^ ^] + ^• 

Set s(a, M,i) = ^ + 2^ + (l + ^)^^'^Q^I^i. Then, s{N , M, ^ = s(/ijv, M, 
^ + 6) and s{N,M,^) = s{uN,M,i + 5'), where 5,5' ^ Q as M,N ^ oo \i 
M-'^{\ogNY/^ converges to as M, iV ^ oo. Thus, (1.27) follows from (4.8) 
and (4.9) and the theorem is proved. D 



5. Applications of the Krawtchouk ensemble 

5.1. Seppdldinen^s first passage percolation model. The Krawtchouk en- 
semble is defined by (1.7) with the weight w{x) = {^)p^q^~^, < x < K, i.e. 
we consider the probability measure 

N 



rKr,N,KAh] = ^-^AMihf TT f f) 



phjqK~hj 
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on {0,...,i^}^, where Z^^k^p = Nll^"^-^^ jj^\k\)^ {pq)^^^-^)l\ The 

first problem where the Krawtchouk ensemble appears is in the simplified first- 
passage percolation model introduced by Seppalainen in [Se2]. Consider the 
lattice N^ and attach a passage time r(e) to each nearest neighbour edge. If 
e is vertical r(e) = tq > 0, and if e is horizontal then T(e) is random with 
P[r(e) = X\=p and P[r(e) = k\= q=l-p, where k > A > 0, < _p < 1. All 
passage times assigned to horizontal edges are independent random variables. 
Hence, all randomness sits in the horizontal edges. The minimal passage tim,e 
from (0, 0) to (k, I) along nearest neighbour paths is defined by 

(5.1) r(/c,/) = min^r(e) 

egp 

where the minimum is over all non-decreasing nearest neighbour paths p from 
(0, 0) to {k, I). The time constant is defined by fj,{x, y) = lim.n—^oo-T{[nx], [ny]). 
(The existence of the limit follows from subadditivity.) In [Se2] it is proved, 
using a certain associated stochastic process, that 



(5.2) fi{x,y) - ^ x__ , __ , / ^^/ /— ,— ^2 



Xx + Toy, if py > qx 

\x + Toy + {k - \){y/qx - y/pyf , if py < qx. 

We will show that the distribution of the random variable T{k, I) relates to the 
distribution of the rightmost charge ("largest eigenvalue") in a Krawtchouk 
ensemble. 

Write M = k,N = /+ 1 and consider an M x A^ matrix W whose elements, 
w{i,j), are independent Bernoulli random variables, P[w{i,j) = 0] = q and 
P[w{i,j) = 1] = p = 1 — q. Let IIm,n be the set of all sequences vr = 
{{k,Jk)}kLi such that 1 < j'l < • • • < Jm < A, i.e. up/right paths in W with 
exactly one element in each row. Introduce the random variable 



(5.3) L{W) = max < 



y~^ w{i,j); vr G IIm.at 



> . 



Write p = 1/q—l, so that q = (l+p) ^ andp = p(l+p) ^. It is straightforward 
to show that 

(5.4) T{k,l) = lTo + kK-{K-X)L{W). 

Proposition 5.1. Let L{W) be defined by (5.3) with W an M x N O-l- 
matrix with independent Bernoulli elements w{i,j), the probability off being p. 
Then, 



(5.5) P[L{W) <n]= Pk 



i,N,N+M-l,p 



max hj < n + N — 1 
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Proof. Interpreting the formula (7.30) in Theorem 7.1 in [BRl] in the 
appropriate case, we get 

e(\) 
(5.6) P[L(T^)<n] = (l + p)-^^ Y. d^iM)dyiN)l[p^^, 

e{X)<n 

where A' is the partition conjugate to A, A^ is the length of the k^^ column 
in A, and dx{M) is the number of semi-standard tableaux of shape A with 
elements in {1, ... , M}; if ^(A) < M, dx{M) is given by (4.2). The proof of 

(5.6) is based on the RSK-correspondence between 0-1 matrices and pairs of 
semistandard Young tableaux (P, Q) where P has shape A and Q has shape 
A', see [Fu], [St]. Set 

^m{N) = {A G f^M ; A^ > Ai > . . . Am > 0}. 

Since dx{M) = if l{\) > M and dyiN) = if Ai > iV, (5.6) can be written 
as 

(5.7) P[L{W) < n] 

=(i+p)-"MnJF n^ E ^M(A)y.(A')n/^ 

l{X)<n 

Note that A G ^m{N) if and only if A' G ^n{M) and l{\) = \[. 
Lemma 5.2. If jx ^ Qn{M), then 

(N+M-l \ N 

Proof. One way to prove (5.8) is to use the fact that VM{tJ'')WM{lJ'') 
= yM(M)^M(/ti) by the hook formula for /^; compare with (1.13) and Lemma 
3.1. We will give another proof. Set Si = Hi + N + l — i, l<i<N and 
Vj = N + j - fi'j, l<j<M. Then, 

(5.9) {si,...,SN}U{ri,...,rM} = {l,---,N + M}. 

To see this, notice that since 1 < Si,rj < A^ + M it suffices to show that 
Si 7^ rj for all i,j. Looking at the //-diagram one sees that fii + fi'j < i + j — 2 
or Hi + n'j > i + j, which implies Sj / rj. 

Let nfc = 1 if /c G {si, . . . , sj\[} and n^ = if A; G {ri, . . . jVm}, k = 
1,...,N + M. Then, by (5.9), 

(5.10) Vm{ix')= n (/-fe)(i-"'=)(i-"'). 

l<fe</<Af+M 
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Now, 

l<k<l<N+M 

N N+M N 

n il-kr = ll n il-Sj) = ll{N + M-s,y. 

l<k<l<N+M j=l l=Sj+l j=l 

and 

N Sj~l N 

l<k<l<N+M jr=l fc=l j=l 

Inserting this into (5.10) gives the formula (5.8). The lemma is proved. D 

We can now finish the proof of the proposition. If we write ^ = A', we see 
from (5.8) that (5.7) can be written as 



(5.11) 



P[LiW) < n] = (1 + pr^'' n' ^:^^ 

i=o ^' 



fil<n 

As usual we introduce the new coordinates hj = fij + N — j. Then, using 
p = 1/q — 1, we obtain 

P[L{W) < nj - ^ J^ j! ((iV + M-1)!)^ 



E A^(^)^n ("""->% 



hi^N+M-l-hj 



max(?ij ) <n+ A'^— 1 

which completes the proof. D 

Using Proposition 5.1 we can prove a limit theorem for the first passage 
time T{k, I). The result should be compared with Remark 1.8 and Conjecture 
1.9 in [Jo3]. 

Theorem 5.3. If p{x,y) is given by (5.2), 
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and py < qx, then 

^^ T{[nx],[ny])-Mx,y) < ^i , _ ^(_^)^ 
n-^oo a{x,y)n'^'-^ 

where F{t) is the Tracy-Widom distribution (1.5). 

Proof. The proof uses (5.4) and Proposition 5.1 and is analogous to the 
proof of Theorem 1.7, the difference being that we now need the analogue 
of Lemma 4.1 for the Krawtchouk polynomials. This can be obtained from 
a steepest descent analysis of the integral formula for these polynomials in 
much the same way as in the analysis of the Meixner polynomials in Section 
5 of [Jo3]; see [Jo4] for some more details. The time constant is related to 
the right endpoint of the support of the equilibrium measure associated with 
the Krawtchouk ensemble, and the constant a{x, y) comes out of the steepest 
descent argument. We can also get large deviation results by using the general 
results of Section 4 in [ Jo3] . D 

5.2. The Aztec diamond. We turn now to the relation between the 
Krawtchouk ensemble and domino tilings of the Aztec diamond introduced 
by Elkies, Kuperberg, Larsen and Propp in [EKLP]. The definitions are taken 
from that paper and the papers [JPS] and [CEP] where more details and pic- 
tures can be found. A domino is a closed 1 x 2 or 2 x 1 rectangle in M^ with 
corners in 7? , and a tiling of a region i? C M^ by dominoes is a set of dominoes 
whose interiors are disjoint and whose union is R. The Aztec diamond, An, of 
order n is the union of all lattice squares [m, m + 1] x [1,1 + 1], m,l € 7j, that 
lie inside the region {{x, y) ; |x| + \y\ < n + 1}. It is proved in [EKLP] that the 
number of possible domino tilings of An equals 2"'"'+^)' ^. Color the Aztec dia- 
mond in a checkerboard fashion so that the leftmost square in each row in the 
top half is white. A horizontal domino is north-going if its leftmost square is 
white, otherwise it is south-going. Similarly, a vertical domino is west-going if 
its upper square is white, otherwise it is east-going. Two dominoes are adjacent 
if they share an edge, and a domino is adjacent to the boundary if it shares an 
edge with the boundary of the Aztec diamond. The north polar region is de- 
fined to be the union of those north-going dominoes that are connected to the 
boundary by a sequence of adjacent north-going dominoes. The south, west 
and east polar regions are defined analogously. In this way a domino tiling 
partitions the Aztec diamond into four polar regions, where we have a regular 
brick wall pattern, and a fifth central region, the temperate zone, where the 
tiling pattern is irregular. 

Consider the diagonal of white squares with opposite corners Q^., k = 
0, . . . , n -|- 1, where Q^. = (— r -\- k,n -\- 1 — k — r), r = 1, . . . ,n. A zig-zag path 
Zr in An from Qq to Q^+i is a path of edges going around these white squares. 
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When going from Q\ to QJc+i ^^ ^^'^ 8° either first one step east and then one 
step south, or first one step south and then one step east. A domino tihng on 
An defines a unique zig-zag path Z^ from Qq to Q^n+i ^^ ^^ require that the 
zig-zag path does not intersect the dominoes. Similarly, we can define zig-zag 
paths from Pq = (— r, n — r) to P^ = (n — r, — r) going around black squares. 
We consider random tilings of the Aztec diamond, where each of the 
2n(n+i)/2 possible tilings have the same probability. This induces a proba- 
bility measure on the zig-zag paths. Consider a zig-zag path in An from Qq 
to Qn^i around white squares. Let hr < • ■ ■ < hi he those k for which we go 
first east and then south when we go from Q^ to Ql.:i, k = 0, . . . ,n; there are 
exactly r such k if the zig-zag path comes from a domino tiling, [EKLP] . Call 
this zig-zag path Zr{h). 

Proposition 5.4. Let {hi, . . . ,hr} C {0, ...,n} be the positions of 
the east /south turns in a zig-zag path Zr{h) in the Aztec diamond An from 
(— r, n + 1 — r) to (n + 1 — r, — r) around white squares. Then, the probability 
for this particular zig-zag path is 

(5.12) P[Z,(/i)]=PKr,.,„,i/2[/i]. 

// {hi, . . . , hr} C {0, . . . , n — 1} are the positions of the south/ east turns in a 
zig-zag path Z'^ih) in An from (— r, n — r) to {n — r, — r) around black squares, 
then 

(5.13) P[Z'M]=^Kr,r,n-lM2[h]. 

Proof. Let Ur{h) be the number of possible domino tilings above Zr{h) in 
the Aztec diamond. From the arguments in [EKLP], see also [PS], it follows 
that 

(5.14) ^Y,(/i) = 2^('-i)/2 n ^^1^^. 

-'■-'■ j — i 

l<i<j<r -^ 

Let fci < • • • < kn+i-r be defined by 

{ki,...,kn+l-r} = {0, ...,n} \{hi,...,hr}. 

If Cr{h) is the number of domino tilings of the region below Zr{h) in An, then, 
using the symmetry of the Aztec diamond, we see that 

(5.15) £,(/i) = 2("+i-'-)("-'-)/2 Yl ^'~^' 



1 — I 

l<i<j<n+l~r ■' 



Thus, the probability for a certain zig-zag path Z,. = Zr{h), specified by h, is 
(5.16) 

9(n+l-r)(n-r)/2+r(r-l)/2 h- h- h—h- 

nz.m- — ^situt^ n 7^ H ^- 

l<«<i<i' l<i<j<n+l-r •' 
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If we let /ij = ^i + r — i, 1 < i < r and kj = r + j — 1 — n'j, 1 < j < n + 1 — r , 
then n and n' are conjugate partitions; compare with the proof of Lemma 5.1 

(iV = r, M = n + 1 — r, Si = hi + 1, Vj = kj + 1). We see from (5.16) that 



P[Zrih)] = 2-("+i-'^> 



'r— 1 ^ \ I n—r -. \ 



where ^ G Q.j.{n + 1 — r). We can now apply Lemma 5.1, which gives 

^[^■■("'i = TSF n — !-^'.(M) n („^. ^ , _ J ^. 

Now, hi = Hi + r — i, so we obtain 

(5.17) P[z,(.),.?i^n(!i^A„,.y^n(,;)^. 

which is the Krawtchouk ensemble. Note that in (5.17) the order of the /ij's 
is unimportant, so we can let {/ii, . . . , hr] CI {0, . . . , n} be the (unordered) 
positions of the east/south turns. A completely analogous argument applies to 
the zig-zag paths in An from Pq to P^ around black squares. This completes 
the proof. D 

It is proved in [JPS] that, with probability tending to 1 as n ^ oo, the 
asymptotic shape of the temperate zone is a circle centered at the origin and 
tangent to the boundary of the Aztec diamond {the arctic circle theorem). 
This can be deduced from Proposition 5.4 and the general results in Section 4 
of [Jo3]. The arctic circle is determined by the endpoints of the support of the 
equilibrium measure (or the points where it saturates). Also, from Theorem 
5.3, we see that the fluctuations of the temperate zone around the arctic circle 
is described by the Tracy- Widom distribution. This can also be deduced from 
the fact, derived in [JPS], that the shape of a polar region is related to the 
shape of a randomly growing Young diagram. The growth model obtained is 
exactly the discrete time growth model studied in [Jo3], and we can apply the 
results of that paper. See [Jo4] for more details. 

Finally, we will shortly discuss another random tiling problem related to 
plane partitions using the combinatorial analysis by Cohn, Larsen and Propp 
in [CLP]. For more details and pictures see the paper [CLP]. Plane partitions 
in an a X 5 X c box can be seen to be in one-to-one correspondence with tilings of 
an a, 6, c-hexagon with unit rhombi with angles 7r/3 and 27r/3, called lozenges. 
An a, 6, c-hexagon has sides of length a,b,c,a,b,c (in clockwise order), equal 
angles and the length of the horizontal sides is b. If the major diagonal of 
the lozenge is vertical we talk about a vertical lozenge. Consider the uniform 
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distribution on the set of all possible tilings of the a, b, c-hexagon with lozenges, 
which corresponds to the uniform distribution on all plane partitions in the 
axbxc box. For simplicity we will now restrict ourselves to the case a = b = c. 
A horizontal line k steps from the top, k = 0, . . . ,a will intersect the vertical 
lozenges at positions hi + 1, . . . ,hk + 1, < hi < • • • < h^ < a + k — 1, 
otherwise it passes through sides of the lozenges. A random tiling induces 
a probability measure on the sequences h = (hi, . . . ,hk). Interpreting the 
formulas in Theorem 2.2 in [CLP] we see that the probability for h is 

where Zk^a is a constant that can be computed explicitly. Note that the 
measure is symmetric in the hi's so we can regard (5.18) as a measure on 
{0, . . . ,a + k — l}*^. Thus, again we get a discrete orthogonal polynomial en- 
semble, this time with the weight 

on {0, . . . , N}, with a = (3 = a — k and N = a + k — 1. The orthogonal 
polynomials for this weight are the Hahn polynomials, [NSU], so (5.18) should 
be called the Hahn ensemble. If we do not have a = b = c we will again get a 
weight function of the form (5.19) but with different values of a, (3 and A^ and 
with a different number of particles. This model is further discussed in [Jo4], 
but to obtain the Tracy- Widom distribution in this model is more complicated 
due to the fact that it is less straightforward to compute the asymptotics of 
the Hahn polynomials. 
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